ervices 


Technical  Information  Agency 


Reproduced  by 

DOCUMENT  SERVICE  CENTER 

KNOTT  BUILDING  DAYTON,  2,  OHIO 


This  document  is  the  property  of  the  United  States 
Government.  It  is  furnished  for  the  duration  of  the  contract  and 
shall  be  returned  when  no,longer  required,  or  upon 
recall  by  ASTIA  to  the  following  address: 

Armed  Services  Technical  Information  Agency,  Document  Service  Center, 
Knott  Building,  Dayton  2,  Ohio. 


NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A  DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


NAVY  DEPARTMENT 
THE  DAVID  W.  TAYLOR  MODEL  BASIN 
Washington  7,  D.C. 

CONTRA-ROTATING  OPTIMUM  PROPELTERS  OPERATING 
IN  A  RADIALLY  NON-UNIFORM  WAKE 

by 

Ho  W*  Lerbs 


F 


May  1955 


Report  941 


CONTRA -ROTATING  OPTIMUM  PROPELLERS  OPERATING 


IN  A  RADIALLY  NON-UNIFORM  WAKE 


by 


H.W.  Lerbs 


May  1955 


Report  941 


ii 

TABLE  OF  CONTENTS 

page 

INTRODUCTION . ,■ .  1 

GENERAL  CONSIDERATIONS .  2 

COMPONENTS  OF  THE  RELATIVE  VELOCITY  AND  OF 
THE  FORCES .  3 

CRITERION  FOR  CONTRA-ROTATING  PROPELLERS  OF 
MINIMUM  LOSS  OF  KINETIC  ENERGY .  7 

DETERMINATION  OF  THE  CIRCULATION  FUNCTION  FOR 
OPTIMUM  FLOW .  10 

CORRECTIONS  FROM  THE  WAKE  DIFFERENCE  AND  FROM  THE 
RACE  CONTRACTION .  12 

DESIGN  PROCEDURE .  17 

OPTIMUM  CIRCULATION  FUNCTION  FOR  UNIFORM  INFLOW 
AND  COMPARISON  WITH  TREODORSEN  » S  RESULTS .  20 

REFERENCES . 22 


iii 


T 

Q 

P 

F* 

z 

D  -  2R 
x  =  r/R 
d 

S' 

vs 

w 

w 

o 

V  =  Vg(l-W) 

t 

t„ 


w 


a 


wL 


n 

CJ 


r  v  /jrnD 


Cip  = 


Ih2TTV3 


NOTATION 


Total  thrust  ' 

Tangential  force 
Total  power  j 


in  non-viscous  flow 


Pitch 


Number  of  blades  of  one  of  the  propellers 
Diameter 

Non-dimensional  radial  coordinate 
Distance  of  the  two  propellers 
Contraction  ratio 
Ship  speed 

Effective  wake  factor 
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Local  velocity  of  advance 

Thrust  deduction  factor 
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Power  loading  coefficient  in  non-viacous  flow 

Lift  coefficient  of  any  section 

Drag-lift  coefficient  of  any  section 

Bound  circulation  of  one  blade 

Non-dimensional  bound  circulation  of  one  of 
the  propellers 

Propulsive  coefficient 

Average  factor 

Distance  factor 

Induction  factor 

A  bar  means  an  average  value  taken  either 
circumferentially  or  radially 


Indices : 
1 
2 
3 

i 
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Front  propeller 
Rear  propeller 
Self  induction 
Interference  induction 
Hub 


INTRODUCTION 


Known  theories  of  a  contra-rotating  propeller  are  either 
restricted  to  uniform  Inflow  (1)#  or  Include  arbitrary  assumptions 
concerning  both  the  applicability  of  the  Goldstein  function  and 
the  orientation  of  the  resultant  induced  velocity  relative  to  the 
free  vortex  sheets  (£,  3).  These  assumptions  are  avoided  in  the 
following  considerations  which  make  use  of  the  so-called  induction 
factor's  of  vortex  sheets o 

A  criterion  for  optimum  flow,  expressed  in  terms  of  the 
direction  of  the  free  vortex  sheets,  is  obtained  from  first  order 
considerations.  This  criterion  leads  to  a  non-linear  integral 
equation  for  the  optimum  circulation  or,  approximately,  to  a  set 
of  non-linear  algebraic  equations  for  the  Fourier  coefficients  of 
the  circulation.  For  uniform  inflow,  the  free  vortex  sheets 
become  of  a  true  helical  shape  and  the  equations  for  the 
circulation  reduce  to  a  linear  system. 

A  design  method,  which  follows  from  the  considerations,  is 
outlined  taking  approximately  into  account  the  effects  arising 
from  the  difference  of  the  wake  at  the  propeller  discs  and  from 
the  contraction  of  the  r.ce  between  them. 

Finally,  the  optimum  circulation  obtained  by  Theodorsen  for 
uniform  Inflow  by  means  of  an  electrical  analogy  (1)  is  compared 
with  the  result  from  the  developed  relations. 


The  numbers  in  parentheses  refer  to  the  list  of  references  on 
page  22 
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GENERAL  CONSIDERATIONS 

The  flow  at  either  propeller  of  a  contra-rotating  pair 
arises  both  from  self  induction  and  from  interference  induction 
(Figure  1). 

At  the  front  propeller,  the  interference  causes  an  increase 
of  the  axial  component  of  the  relative  velocity.  There  is  no 
interference  effect  on  the  tangential  component,  as  follows  from 
Stokes'  law. 

At  the  rear  propeller,  both  the  axial  and  the  tangential 
component  of  the  relative  velocity  are  increased  by  interference. 

The  following  assumptions  are  made  when  establishing 
expressions  for  the  flow  components  and  for  the  generated  forces: 

(1)  The  shape  of  the  vortex  sheets  is  determined  by  the 
relative  flow  at  the  disc.  This  assumption  requires  both  that 
the  propellers  are  moderately  loaded  and  that  the  wake  does  not 
appreciably  change  in  the  axial  direction. 

The  first  holds  sufficiently  if  secoi  1  and  higher  powers 
of  the  induced  velocity  may  be  neglected,  tlu  second  if  the 
frictional  part  of  the  wake  is  predominant. 

(2)  The  self  and  mutual  interference  of  he  free  vortex 
sheets  may  be  neglected. 

(3)  The  blades  may  be  replaced  by  lifting,  lines  and 
corrections  for  lifting  surface  effects  may  be  introduced 
afterwards. 

(4)  The  induced  velocity  components  may  be  considered 
time  averages,  i.e.,  oscillations  in  both  angle  of  attacL  and 
magnitude  of  velocity  are  not  taken  into  account. 

(5)  The  bound  circulation  may  be  represented  by  a  Fourier 
sine  series. 

(6)  The  propellers  ha' c  an  equal  number  of  blades.  This 
latter  assumption  can  be  easi’y  removed  if  deemed  necessary. 
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COMPONENTS  07  TITS  RELATIVE  VELOCITY  AND  OF  THE  FORCES 

From  Figure  1,  the  components  of  the  relative  velocity  at 
a  blade  element  are  for  the  front  screw; 

axial  =  v  +  (wag)1  (w&i)1 

tangential  =  tfr  -  (wtg) 


and  for  the  rear  screw; 

axial  =  v  +  (w  )„+  (w  . )„ 
as  2  ai'2 

tangential  =  A»r  -  j\'wfca)g  -  (wtl)2] 

The  interference  velocities  are  considered  circumferential  averages, 
from  assumption  4.  These  averages  are  obtained  when  multiplying 
the  maximum  value  of  the  respective  induced  velocity  component, 
which  is  induced  at  a  bound  vortex,  by  its  average  factor  f  and 
when  introducing  a  second  factor  g  to  represent  the  change  of  the 
average  in  the  vicinity  of  the  disc.  Correspondingly,  the 
interference  velocities  may  be  expressed  by  the  self-induced 

(ea)2l=  (%T}2  t  '  (sa's] 


(Sa>i]r  (was>!  [l  ♦  (Sa)J 

f  (et)i]a  2  [l  ♦  (6t)J 

The  factor  2  within  the  last  relation  arises  from  the  discontinuity 
of  the  tangential  component  of  the  self  induced  velocity. 

To  approximately  determine  the  factors  ga  each  propeller  is  replaced 
by  a  uniformly  loaded  sink  disc  and 

( Sa  ^ i  t  ( Sa )  2 


velocities  as  follows; 

(w  .)  =  <w  )„  (f  )„  f1  - 

ai  l  as  2  a  2  L 

(*ti>i  =  ° 

(wai>2  =  (was^i  (fa>l  \}  * 
(wti)2  =  2  (ft)i  [l 


is  assumed,  i.e.,  symmetry  of  the  gradient  of  the  axially  induced 
flow  relative  to  the  disc  is  assumed.  From  the  known  potential 
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of  a  sLnk  dLsc  it  develops  that 


1  „  ,  5  1  .  35  1  „ 

S4%  *  w  S6p«  Tii  S8  p6  ' 


6] 


rfl  1  3 

ga  =  1  -  cos  ©^  b2-e 

*  sm  e  |  ^2-^  56  p4  ^  58  Pe  -  ♦  •••] 


if  b  > 1  and 


ea  :  1  -  C03  ej^Px  -  b  P2  t  |  b3  P4  -  |  b5  P6;*  -  ...j 
f  sin  9  |px  -  I  b  pj  t  i  b3  P;  -  i-b3  Pi  ♦  - 

if  b<  1. 

The  functions  Pn  are  the  Legendre  polynomials  of  the  argument 


cos  © 


=  d/R 


where 


b  =  fa2  -  (d/R)2 


and  where  d  represents  the  axial  distance  of  the  two  propellers. 

The  prime  means  the  derivative  of  P  relative  to  9.  The  factor 
g  is  represented  on  Figure  2  as  a  function  of  x  and  d/R«  The 
diagram  indicates  that  the  change  of  the  axial  component  is  rather 
great  in  the  vicinity  of  the  disc.  This  behaviour  may  necessitate 
replacing  the  approximations  [l}  by  more  rigorous  expressions 
taking  into  account  the  non-uniiorm  loading  of  the  disc. 

The  factor  g.  depends  on  contraction  and  is  assumed  to  be  zero  for 
the  present.  This  assumption  will  be  corrected  on  page/3,, 
CORRECTIONS  FROM  THE  WAKE  DIFFERENCE  AND  FROM  THE  RACE  CONTRACTION. 

The  average  factor  ft  follows  from  Stokes'  law.  If  (zT  ) 
represents  the  total  bound  circulation  of  one  of  the  propellers. 
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this  law  yields  the  result 

2  rir(2wtg)  ft  =  2  r7T(5wts)  =  z  T 

from  v/hich  it  follows  that 

[2]  ft  ={hi}  S  ZP 

wts  4r?7(wts} 


The  average  factor  of  the  axial  component,  viz.. 


is  not  yet  reliably  known.  Work  for  its  determination  has 
been  initiated  on  a  basis  of  the  Biot-Savart  integral  for  semi- 
infinite  helical  vortex  lines.  ?’e  put  tentatively  until  accurate 
results  are  available: 


£0  fa  !  ft 


Prom  these  relations,  it  follows  that  tu<.  pitches  of  the  free 
vortex  sheets  are  slightly  different  for  *he  front  and  rear 
propeller.  To  simplify,  vortex  sheets  of  equal  shape  are  assumed 
the  pitch  of  which  may  be  taken  as  the  mean  value  of  )t  and 

(tan^1)2. 


A  further  simplification  arises  from  the  requirement  that  the 
torque  of  the  front  propeller  equals  that  of  the  rear  propeller. 
Prom  the  law  of  moment  of  momentum  this  condition  is  satisfied 
if 


[4]  (zr)1  -  (zf  )2  «  ur ) 

The  assumption  involved,  viz.,  that  the  propellers  have  equal 
diameters  will  be  removed  in  CORRECTIONS  PROM  THE  WAKE  DIFFERENCE 
AND  PROM  THE  RACE  CONTRACTION  on  page  12  0  The  bound  circulations 
of  the  two  propellers  being  equal  functions  of  x,  the  circulation: 
of  the  free  vortex  sheets  are  also  equal.  Since  the  slight 
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difference  between  the  pitch  of  these  sheets  is  disregarded, 
the  respective  components  of  the  induced  velocities  become  equal. 
That  is 


Ksh  -  (was ^2  =  Wa 
(wts)l  s  (wts}2  =  Wt 


With  these  simplifications,  the  following  expressions  for  the 
components  of  the  relative  velocity  are  obtained  ~ 


/  fcr  the  front  screw; 


axial  3  v  +  w. 


1  +  fa  (1 


->] 


M 


tangential  ■  6)r  -  w^. 

and  for  the  rear  screw: 


axial  =  v  +  w. 
tangential  3  «r  -  ’ 


[l  *  fa  (1  ♦  ga)] 
■  (  1  “  2  ft) 


Introducing  a  non-dimensional 

[«]  o. — sE — 

L  7TD  v3 


circulation,  viz., 


the  thrust  coefficient  and  the  power  input  coefficient  of  the 
contra-rotating  propeller  are  obtained  from  the  law  of  Kutta- 
Joukowsky  as  follows; 


^  (1  -  ft) 

Vs 


dx 


c 


a 

P 


p  « 

col 

(1  -  w)  t  !^(1  t  fa)  xdx 

As  j 

xh 

vs 

w 
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The  loading  coefficients  so  ascertained  are  slightly  different 
from  the  given  values  Cj  or  cp,  respectively,  from  reasons  which 
are  discussed  in  CORRECTIONS  PROM  THE  WAKE  DIFFERENCE  AND  FROM  THE 
RACE  CONTRACTION  on  page  12.  Within  these  expressions,  the  average 
factors  are  determined  from  the  relation 

[9]  ft  •  - 2 - -  ;  fa 

2x  wt 
Vs 


Further,  the  self  induced  velocities  wa  and  wt  are  related  to  G  by 
the  following  expressions; 


( 


1 

2  z 


(  ia 

l  dxQ  x-x0 


dx0 


=  1 
z(l-xh) 


m 


r*  ha 
Gm  hm 


1 

2z 


dG  1 
dx£  x-x0 


1t 


dx„ 


7  “77 - T  m=l 

z(l-xh) 


m 


Gm 


The  functions  ia  and  i  are  the  induction  factors  of  the  free 
vortex  sheets  as  definid  in  (4).  The  functions  enter  when  the 
integral  on  the  right  hand  side,  which  is  an  improper  integral,  is 
evaluated.  These  functions  are  related  to  the  induction  factors, 
the  relations  having  been  developed  in  (4).  .  By  the  factors  G*  , 
the  coefficients  of  G  are  denoted  when  G  is  represented  by  a 
Fourier  sine  series  within  the  interval  x  9  xj{  and  x  ■  1. 


CRITERION  FOR  CONTRA-ROTATING  PROPELLERS  OF  MINIMUM  LOSS  OF 
KINETIC  ENERGY 

The  e qua  t  i  ons  £7}  t o  [lo]  together  with  an  assumption  on  the 
thrust  deduction  factor  suffice  to  treat  the  variational  problem 
of  a  contra-rotating  propeller  operating  in  a  given  wake  field, 
viz.,  to  determine  the  coefficients  of  the  Fourier  series  for 
G(x)  such  that  the  useful  work  becomes  a  maximum  value  for  a  given 
power  input.  However,  this  problem  of  determining  the  optimum 
circulation  function  in  a  direct  way  is  laborious  for  a  finite 
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number  of  blades  and  Is  not  solved  he»e.  Instead,  a  relation 
between^,  and  y#  is  deduced  for  optimum  flow  from  which  the 
circulation  function  follows  afterwards. 

From  considerations  on  moderately  loaded  single  propellers, 
it  is  known  that  the  form  of  the  optimum  relation  between and 
^is  independent  of  the  number  of  blades.  This  enables  us 
to  consider  propellers  with  infinitely  many  blades  by  which  the 
deduction  of  the  relation  is  greatly  simplified.  We  apply  the 
rule  that  the  elementary  propulsive  coefficient  which  is  related 
to  an  increment  of  circulation  is  independent  of  the  radius  for 
optimum  flow,  see  (4).  If  the  increments  of  the  forces  which 
arise  from  the  increment  of  circulation,  Aun  are  denoted  by  Afd'-i) 
and  A(dQ),  respectively,  the  rule  reads; 

A(dT)  v_  1  -  t(x)  =  k2  , 

A(C'Q)  *>r  1  -  w(x) 

the  constant  k  being  independent  of  the  radius  x. 

To  establish  A(dT)  and  A(d$)  we  write  from  the  law  of  Kutta- 
Joukowsky  for  infinitely  many  blades,  i.e.,  for  fa  s  f^  =  Is 

dT  ■  dT]_  t  dT2  =  2(zp)  fuv  sr 
dQ  -  dt^  +  dQ2  =  2(zf)  />  (v  t  2  w&)  dr 
It  follows  that 

A(dT)  =  2j>av  A(zD  <ir 

and  that 

A(dQ)  =  2J>  f  (v  ♦  2  wa)  +  (zf'2  ,_Awa - 1  A(zT)  dr 

L  A(zT)  J 

A  relation  between  wa  and  (z P)  is  obtained  from  the  law  of  momentum. 
Neglecting  terms  of  second  power  of  wa  this  relation  reads; 

w  &  1  "  t ) 

a  v  4r  Jo 
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Then, 


A(dQ)  -  2/(v  «■  4  wa)  A(zt1)  dr 
and  the  optimum  condition  becomes 


tor  v  =K2  1  -  w(x) 

v  +  4wa  6Jr  1  -  t(x) 


Since  the  axial  velocity  in  the  slipstream  far  aft  equals 


(1  -  w)  +  4  w 


the  condition  means  that 


t  (x 


IL  I"1 

K2'  L 


>] 


If  the  change  of  t(x)  over  the  radius  is  considered  small  so  that 
t(x)  4  t0,  the  condition  requires  that  the  axial  velocity  within 
the  slipstream  far  aft  is  independent  of  the  radius  for  infinitely 
many  blades.  These  relations  are  approximations  which  include 
assumption  1,  viz.,  that  the  change  of  the  wake  in  the  axial 
direction  may  be  neglected. 


To  express  the  minimum  condition  by  the  pitch  of  the  free  vortex 
sheets  we  consider  the  average  of  the  pitches  of  the  two  sheets. 
To  the  first  order,  this  average  for  infinitely  many  blades 
amounts  to 


5  v  ♦  g  »a 

tor 

Since,  to  the  first  order 

v 

tor 


tor 


v  ♦  4  wn 


0)  r 


v  ♦  2 


y 


the  opcimum  condition  may  be  written  as 


[ll]  tan,*? 


tan/^ 
"Tc 


7lh 

/  1  -  w 


-  t(x) 


U) 


p(x) 

k 
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That  Is,  the  mean  pitch  of  the  vortex  sheets  of  a  contra¬ 
rotating  propeller  follows  to  the  first  order  the  same  law  lr 
the  optimum  case  as  established  previously  for  a  single  wake 
adapted  optimum  propeller,  see  (5).  However,  the  optimum 
circulation  function  of  the  contra-rotating  propeller  differs 
from  that  of  a  single  wake  adapted  propeller  as  seen  from  the 
following  section<> 

DETERMINATION  OF  THE  CIRCULATION  FUNCTION  FOR  OPTIMUM  FLOW 

Analogous  to  single  propellers,  the  optimum  relation  as 
expressed  by  fill  is  considered  independent  of  the  number  of 
blades o  Then,  the  optimum  circulation  function  G(x)  is  obtained 
for  a  finite  number  of  blades  when  equating  two  expressions  for 
tan^  viz.,  the  minimum  condition  [llj  and  the  geometric  relation 

v 

tan^  §  v  ♦  wa  (1  »  fa)  -  «r  wt(l  -  ft) 
oiv  -  2  wfc  (1  -  ft) 


This  latter  expression  is  a  first  order  approximation  for  the 
mean  value  of  (tan/i  )i  and  (tany/.Jp  for  a  finite  number  of  blades 
which  follows  from  [5J. 

Introducing  0]  for  the  average  factors  and  equating,  the  relation 

(1  -  w)  -  ±_  P(*)l  +  £_  [2a  f  (1  -  w)  -  2  p(x)  ~1 

^s  k  J  2x  Lw^  x  kj 

=  zt  (1  _w)-  2  PW.1-  22 

vs  L  x  k  J  vs 


is  obtained  within  which 

p(x)  =  tany^fUl HEls  J(  1-w )  ( 1-tT 

/  \1  -  w(x)  X  " 

Is  considered  a  known  function  of  the  radius. 

Expressing  w&  and  w.  by  the  Integral  representation  of  [10]  an 
integro-dif f erential  equation  is  deduced  for  G.  An  approximate 
solution  of  this  equation  becomes  possible  when  utilizing  the  sums 
of  [1(3  for  wa  and  w^.  This  leads  to ‘.the  following  system  of 
algebraic  equations  for  the  Fourier  coefficients  G*  of  G: 
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This  relation  is  satisfied  at  m  stations  x  in  order  to  obtain 
m  equations  for  m  coefficients  of  the  Fourier  series! 

G  =  G*  sin  m  <> 
m-1  m 


where 

x  -■  7(1  *  xh)-  7  (1  -  xh)  cos  f 

On  the  right  hand  side  of  [12],  the  ratio  wa /wt  is  not 
expressed  by  the  sums  of  equation  [lO^ .  The  reason  is  to  avoid 
the  numerical  __  solut  ion  of  non-linear  equations  for  the  . 

coefficients  G^»  This  can  be  done  when  introducing  for  wa/w^. 
a3  a  first  approximation  the  relation  ' 


wt  tan^  P  ( x ) 


''  «  , 

i.e.,  assuming  the  condition  of  normality  to  be  satisfied. 
This  approximation  may  be  corrected  by  successive  solutions 
of  the  system  [IS]  . 


It  should  be  mentioned  that  the  essential  difference  between  a 
contra-rotating  and  a  single  wake-adapted  optimum  propeller, 
which  is  treated  in  (4),  lies  in  the  equations  for  G  .  These 
equations  are  non-linear  for  the  former  and  linear  f$r  the  latter. 
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CORRECTIONS  FROM  THE  WAKE  DIFFERENCE  AND  FROM  THE  RACE  CONTRACTION 


So  far,  the  finite  distance  of  the  propellers  has  been  taken 
into  account  by  the  distance  factor  g.  The  possible  difference 
between  the  wakes  at  the  positions  of  the  propellers  has  been 
neglected  and  the  radius  of  the  streamlines  has  been  assumed  to  be 
constant .  Both  the  wake  difference  and  the  finite  distance 
necessitate  satisfying  the  equation  of  continuity  within  the  space 
between  the  propellers.  Corrections  on  the  induced  velocity 
components  ensue  and  the  diameters  can  no  longer  be  assumed  equal. 
With  different  diameters,  slight  changes  of  the  bound  circulation 
become  necessary  in  order  to  maintain  the  condition  of  an  equal 
torque  for  each  propeller.  We  introduce  the  wake  factor  wj  in  the 
plane  of  the  front  propeller  and  the  factor  Wg  in  that  of  the 
rear  propeller  so  that  ' 

vx  ~  v3  (1  -  wx) 

V2  =  VS  (1  -  W2) 


We  assume  that  the  wake  does  not  change  on  the  front  propeller, l.e., 
that  the  wake  factor  is  identical  with  the  quantity  w  used  in  the 
i  foregoing  considerations,  Then,  the  results  for  both  the  self 

I  .  induced  velocities  from  [10]  and  for  G  from  [l2]  hold  for  the  front 

propelle.  since  it  will  be  seen  that  the  slight  change  of  the 
i  pitch  of  the  free  vortex  sheets  which  ensues  on  the  front  propeller 

j  from  the  changes  on  the  rear  propeller  may  be  neglected. 

!  Correspondingly,  the  index  1  is  attached  to  these  quantities  and 

!  '  the  non-dimensional  circulation  G-^  is  defined  by 

* 

(  [l4]  G1  =  (zT  }1 _  »  G 

I  TTD-l  Vs 


The  equation  of  continuity  requires  the  race  to  contract,  A 
streamline  which  passes  through  equal  radii  on  both  propellers 
when  ignoring  contraction  will  now  pass  through  the  radius  r2  at 
the  rear  propeller.  We  write 

[15]  Tg  s  rx  ( 1  -  S ) 

and  consider  both  «f  and  d  cT  small  quantities.  The  so  related 

dr^  •  *' 
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radii  are  denoted  "corresponding"  radii. 

The  question  arises  how  both  the  self  induced  velocities  and  the 
interference  velocities,  as  determined  from  [lO]  and[5]  on  a  basis 
of  an  equal  wake  w^,  change  when  the  wake  at  the  rear  propeller  is 
different  from  w^.  The  significant  changes  occur  at  the  rear 
propeller.  Relative  to  its  self  induction  this  problem  is 
essentially  that  of  a  change  of  the  induced  velocity  components 
when  propellers  of  slightly  different  values  of  diameter,  bound 
circulation  and  pitch  of  the  free  vortex  sheets  are  considered. 

Let  the  bound  circulation  be  U?)l  and  (zT)p  on  corresponding  radii 
of  the  front  and  rear  propeller,  respectively.  We  write 


[ie]  (zF)2  =  (zDq  (1  +*) 


Bothtfand  d  <  are  small  quantities.  From  the  integral  representation 

r  -l  drl 

ofUQl  it  follows  for  the  axially  self  induced  component  at  the  rear 
propeller,  e.g.,  that 


if  $  (r^)  in  Ql5J  may  be  replaced  by  its  mean  value  <T .  Then 
aA«\  =  0.  Further,  the  induction  factors  do  not  change 
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appreciably  when  the  pitch  of  the  fre -  vortex  sheets  changes 
slightly  so  that  the  ter-  containing  *Ua)l  may  be  neglected. 


Ir,  addition,  the  product  dty’ ( zT ) 

dr]_ 

both  d  tf  and  (z  T)^  are  small. 


i 

is  small  of  second  order  since 
the  latter  because  of  moderate 


loading.  If^fis  replaced  by  its  mean  valued  the  self  induced 
velocities  at  the  radius  V2  of  the  rear  screw  are  approximately 
related  to  those  at  the  corresponding  radius  of  the  front  screw  by 


(was  ^2  ?  U  ♦/)  U  +<)  (was)l 


(wts)2  f  (1  ♦  /)  U  \2)  (wts)1 


Prom  these  equations,  the  afore  mentioned  slight  change  of  the 
pitch  of  the  free  vortex  sheets  on  the  front  propeller  follows. 

The  reason  for  the  change  is  that  the  interference  from  the  roar 
propeller  is  different  when  the  self  induced  velocities  of  the 
latter  are  different.  However,  this  effect  may  be  neglected 
because  of  the  negligible  variation  of  the  induction  factors 
which  is  involved  so  that  the  result  from  (TlQl  obtained  with 
w^  =  w  Is  sufficiently  accurate  for  the  self  induction  on  the  front 
propeller. 

To  deterr ine  the  interference  velocities  from  the  expressions  for 
the  self  induced  velocities  the  changes  of  th-  factors  f  and  g  must 
be  known.  The  change  of  the  average  factor  f  on  the  rear  propeller 
arises  from  the  changes  of  the  circulation,  of  the  radius  and  of  the 
self  induced  velocity.  It  follows  from  [2j  ,  however,  that  the 
resultant  change  of  f  Is  small  of  second  order  so  chat 


ft  (r2)  =  ft  (rx) 

fa  (r2>  §  fa  <rl) 

The  change  of  the  distance  factors  go  and  g  arises  from  the 
contraction  of  the  slipstream.  For  the  factor  ga,  the  effect  is 
included  in  equations  03  since  these  equations  represent  the 
rigorous  solution  for  the  inflow  to  the  sink  disc.  For  the  front 
screw,  ga  is  determined  at  the  radius  r^  and  for  the  rear  screw  at 
the  radius  r g.  However,  the  accuracy  of  these  equations  is  less  for 
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the  flow  behind  a  propeller  than  ir..  front  because  of  the  assumed 
symmetry  of  the  gradient  of  the  axially  induced  flow.  Because  of 
this,  the  approximation 

ga(r2)  I  ga(rx) 

is  considered  to  be  in  accora  with  the  applicability  of  equations  [lj. 

The  effect  of  contraction  on  the  factor  g  is  established  from 
the  law  that  the  circulation  remains  unaltered  Bn  a  contracting 
streamline . 

The  components  of  the  interference  induction  are  then  related 
to  the  self  induction  as  follows? 


(wai>i 

=  (w — ) 

as  2 

(1 

6a) 

(wti>l 

=  0 

(w  ) 

=  (W~l 

(1 

t  gfl) 

ai  2 

as  i 

a. 

<wti)2 

=  2(w^i 

(1  f  S ) 

Prom  these  expressions  for  the  self  and  interference  induction  the 
components  of  the  relative  velocity  follow  at  corresponding  radii. 
One  obtains 


'for  the  front  screw  (radius  r^) 

axial  *  v1  f  (wa3)1  *  (wal)1 


N 


-  v3  (1  -w1)  +  (was)1  £l  ♦  fa(l  *S)(1*7)  (1  -  ga) 
=  vs  A1 


tangential  -  O  r^  -  (w^.g)^ 
for  the  rear  screw  (radius  rg  corresponding  to  r^) 
\  axlal  “  v2  +  (was^  +  (wai}2 
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Vgd  “  Wa  )  +  (Wag)1  |  (1  »-^5  (1+/)  tfa  (1 


♦  Sa>] 


=  V3  A2 


tangential  ~  O  r2  -  (wtg)2  +  (wtl)2 

=  Olr1  (1  -S)  -  (wt3)1[(l  +4)  (1  +1)  -  2  ft(l  +i) 


] 


where  (wa3)^  an^  (w^)  are  re3U^  from  £lo]  with  w  “  and  r  =  r 

The  contraction  S  (r, )  is  ascertained  from  the  equation  of 
continuity  which  reads  for  an  annular  element  as  follows: 


[' 


1  +  {was  h  + 


J  rldrl  S[V2  +  (was>2  *  (wai>2]  r2dr2 


Introducing  the  expression  £l5^  and  the  relations  for  (w&g)  and  (w^) 

a  first  order  linear  differential  equation  for  &  is  obtained  when 
neglecting  terms  which  are  small  of  second  order-  Determining  the 
constant  of  Integration  such  that  / -  0  for  x-^  -  (x,  )  the  solution 
is  represented  by  the  following  integral:  l 


(wrV  fa  ga[2  M l  \  <)]  dXi 

(l-w2)  +  (wag)1  fa  2  +  (/  +  ga)J 


It  remains  to  establish  a  relation  for  the  function  4?  (r-^)  which 

follows  from  the  condition  of  an  equal  torque  for  each  propeller. 
From  the  law  of  Kutta-Joukowsky  the  balance  requires  that 


1? 


This  equation  is  In  accord  with  the  law  of  moment  of  momentum  which 
necessitates  that  4-  0  if  the  propellers  are  very  close  together  , 
i.e.,  if  both  $  -  0  and  ga  ~  0  - 


The  functions  S  and  <£&re  interrelated  within  [l8]  and[l9],  These 
equations„may  be  solved  in  successive  steps  assuming  in  the  first 
step  f,  S  and  ^  to  be  zero.  The  solutions  for  &  (x^)  tind  (x^)  so 

obtained  lead  to  first  approximations  for  <T  and  <  with  which 
improved  values  ofi^x^)  and  xn  }  follow. 


With  $  being  known  the  change  of  the  diameter  may  be  determined  which 
is  required  by  the  equation  of  continuity.  This  equation  is  satisfied 
if 


( 1  -  /  )x^  „  ^ 


Finally,  the  alteration  of  the  total  thrust  may  be  ascertained  which 
ensues  from  the  changes  under  consideration.  Between  the  thru3t 
coefficient  c*‘  as  follows  from  [7]  and  the  coefficient  crp  obtained 

from  the  law  of  Kut ta- Joukowsky  when  introducing  equations  [l5j  to 
[17]  the  approximation 

H  CT  T  CJ  [l  -  ] 

holds. 

DESIGN  PROCEDURE 

From  the  developed  relations,  a  design  may  proceed  in  the 
following  way  considering  the  total  thrust  a  given  quantity  which 
is  denoted  c  .  The  first  step  is  to  solve  the  algebraic  equations 

[12]  for  the  coefficients  G*  assuming  several  values  of  the  constant 
k  and  to  determine  k  such  tioat  the  thrust  coefficient  c$  is  obtained 
from  [7]  .  The  quantity^*  in_  [7]  is  related  to  the  required  quantity 
cT  by  [»]•  Since  both  <T  and  -£are  unknown  when  beginning  the 
calculation  an  estimate  is  necessary  to  determine  cS.  A  slight 
increase  of  Crp  will  usually  be  sufficient  to  satisfy  [a]. 
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Tne  constant  k  is  approximately  ascertained  by  the 
following  relation,  see  (5): 

k  *  ( PC )  V  H 


so  that 


y^i = 


P(*i  } 


1  ■-  tf 


w(xx) 


These  relations  make  it  possible  to  approximately  determine  k 
and  tan  from  estimated  values  of  both^the  propulsion 
coefficient  end  t0.  The  assumption  t(x)  ~  t  made  in  the 
expression  for  tan^^  can  be  improved  utilizing  results  of  (6). 

However,  the  influence  of  this  assumption  on  the  coefficients  G* 
is  small. 

With  this  fir3t  approximation  for  tany^  the  functions  i  and 

h  are  established  In  a  way  which  is  analogous  to  that  shown  in  (4) 
fSr  a  single  propeller.  Next,  the  equations  [l2j  are  solved  and 
the  circulation  function  G(x),  which  is  identical  with  Gi(x1),  is 
calculated  from  Il3l  •  The  components  of  the  self  induced  velocity 

..  ~  fioj.  -  - - 


follow  from  the  sums  of 


_ _  _  _  __  v--.-  The  average  factors  are  then  known 

from  [9]  and  the  thrust  loading  coefficient  from  [7] .  If  the  K 
difference  between  the  required  and  the  calculated  value  for  c™  is 
appreciable  the  calculations  are  repeated  using  a  slightly  different 
value  for  k<. 

After  sufficient  agreement  between  the  two  values  for  cT  is 
obtained  both  the contraction S and  the  function^ are  found 
from  [18]  and  [l9j  .  The  diameter  D  follows  from  [20J  ,  D]_  being 
assumed  from  the  beginning.  The  pitches  and  the  products  (c *«£)  are 
then  ascertained  from  the  following  relations  which  hold  at  the 


then  ascertained  from  the  following  relations  ’ 
radii  xx-  ry/Ri  and  x2  =  rg/Rg  ”  xx 


respectively: 


n/i]r 


(tany^). 


+^was^i 


(wai)2 


(wts)- 


t 
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and 
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using  incices  1  or  2  on  the  terms  in  brackets  to  obtain  Bj_  or  Bg 
respectively o  ~  * 


vise 


Pvis< 


|v3 
2 


s  R1tr 


lit'1*'1'  2  V  ^  ] 


( sin jf\.)  t1  *  £  > 

tan^ 


(xdx) 


The  former  equation  may  be  used  to  check  whether  the  rumerical 
value  for  the  ratio  ct/(ct)v1sc  ree-.onably  assumed  when  starting 
t.  s  computations .  The  latter  relation  gives  an  estimate  of  the  power 
input  and,  considering  the  integrals  separately,  can  also  be  used  to 
check  whether  the  balance  of  torque  is  maintained  in  viscous  flow. 

As  in  the  design  of  single  propellers,  corrections  on  the  sections 
for  lifting  surface  effect^  (additional  camber  and  additional  angle 
of  attack)  and  for  viscous  flow  (additional  angle  of  attack)  are 
introduced  when  determining  pitch  and  camber.  Also,  the  effect  of 
the  hub  on  the  induced  velocity  may  be  estimated  as  outlined  in  (4). 

OPTIMUM  CIRCULATION  FUNCTION  FOR  UNIFORM  INFLOW  AND  COMPARISON 
WITH  THEODORSEN's  RESULTS 

For  uniform  inflow,  i.e.,  for 

w  a  o,  t  *  0,  p(x)  a 

x 


the  minimum  condition 


tan 


I 


1 

k 


l 


reduces  to 


That  is,  the  mean  pitch  of  the  vertex  sheets  forms  a  true  helical 
surface.  Consequents ,  the  condition  of  normality  is  satisfied  so 
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that  the  relation 


holds  in  this  case.  Then,  the  equations  [l2]  for  the  optimum 
circulation  are  simplified  to  the  following  linear  system  assuming 
xh=  0 ; 


L  = 


<£  (m. 


m 


m.x 


N 

m,x 


) 


where 

L  =  1  •“  A 

Ai 


H 


m,x 


sin(m  <$ ) 
2x 


A 


X 


This  system  is  solved  numerically  for  the  arbitrary  values 
z  -  4,AS  0,80,  =  0.94,  The  result  obtained  is  compared  with 

that  determined  by  Theodorsen  for  hubless  propellers  by  means  of 
an  electrical  analogy  which  is  based  on  the  fact  that  a  velocity 
potential  may  be  reproduced  by  an  electrical  potential  provided 
that,  the  boundary  conditions  are  identical  (1), 

The  following  solution  of  the  set  of  equations  is  obtained} 

gJ  -  ,11907  G4  =,00024 

G*  =  .00844  at  ~  ,00005 

2  5 

G|  =  ,00417 


The'  ensuing  circulation  distribution  G  is  represented  on 
Figure  3,  As  compared  to  the  circulation  following  from 
Theodorsen' s  work,  differences  arise  near  x*0  where  G  goes  to  zero 
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whereas  Theodorsen's  function  remains  finite.  This  difference 
arises  from  the  representation  of  G  by  a  Fourier  sine  series  which 
is  not  suited  to  represent  the  optimum  circulation  of  a  hubless 
contra-rotating  propeller  close  to  the  axis.  However,  for  the 
essential  part  of  the  blade,  between  x  about  0.3  and  1,  the 
agreement  is  considered  satisfactory.  Thl3  result  indicates  that 
the  optimum  circulation  of  a  contra-rotating  propeller  with  a 
finite  hub  is  sufficiently  represented  by  the  series  chosen  for  G 
also  close  to  the  hub. 
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